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ABSTRACT 



(— i Context. Extrasolar giant planets are found to orbit their host stars with a broad range of semi-major axes 0.02 < a < 6 AU. Current theories 
O suggest that giant planets orbiting at distances between =: 0.02 - 2 AU probably formed at larger distances and migrated to their current 
I locations via type II migration, disturbing any inner system of forming terrestrial planets along the way. Migration probably halts because of 
O ' fortuitously-timed gas disk dispersal. 

Aims. The aim of this paper is to examine the effect of giant planet migration on the formation of inner terrestrial planet systems. We consider 
situations in which the giant planet halts migration at semi-major axes in the range 0.13 - 1.7 AU due to gas disk dispersal, and examine the 
effect of including or neglecting type I migration forces on the forming terrestrial system. 

Methods. We employ an N-body code that is linked to a viscous gas disk algorithm capable of simulating: gas loss via accretion onto the central 
star and photoevaporation; gap formation by the giant planet; type II migration of the giant; optional type I migration of protoplanets; gas drag 
on planetesimals. 

Results. Most of the inner system planetary building blocks survive the passage of the giant planet, either by being shepherded inward or 
scattered into exterior orbits. Systems of one or more hot-Earths are predicted to form and remain interior to the giant planet, especially if 
type II migration has been limited, or where type I migration has affected protoplanetary dynamics. Habitable planets in low eccentricity warm- 
lupiter systems appear possible if the giant planet makes a limited incursion into the outer regions of the habitable zone (HZ), or traverses its 
entire width and ceases migrating at a radial distance of less than half that of the HZ's inner edge. 

Conclusions. Type II migration does not prevent terrestrial planet formation. There exists a wide variety of planetary system architectures that 
can potentially host habitable planets. 
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^ 1 Key words, planets and satellites: formation - methods: N-body simulations - astrobiology 

\ 1 . Introduction. (e.g. Lin & Ida 1997; Ford et al. 2001; Papaloizou & Terquem 

2001; Marzari & Weidenschilling 2002), a combination of mi- 
Giant planets are thought to form in the cool, outer, regions of gration mi scatter i ng (Adams & Laughlin 2003; Moorhead & 
a protoplanetary disk (e.g. Pollack et al. 1996; Papaloizou & Adams 2 005), or migration along with eccentricity excitation 
Nelson 2005; Boss 2000), in roughly the region where Jupiter from the disk ( Pa p a loizou et al. 2001 ; Goldreich & Sari 2003; 
and Saturn are found in our solar system. However, numerous Qgilvie & Lubow 2003; Moorhead & Adams 2008). 
giant exoplanets have been found orbiting solar-type stars well 

inside the approximate position of their nebular snowline with In the case of migrating planets, the mechanism that ter- 

semi-major axes from ~ 3 AU down to just a few stellar radii minates the migration and strands exoplanets at their present 

(Butler et al. 2006). The most extreme examples of these are the orbital radii is unknown. Migration-halting mechanisms that 

so-called 'hot-Jupiters', orbiting within 0.1 AU and account- might work when the planet ventures close to the central star 

ing for about a quarter of the known giant exoplanet inventory, include tidally-induced recession caused by the star's rotation 

Planetary migration may provide the best explanation for the or Roche lobe overflow and mass loss to the star (Trilling et al. 

presence of the hot-Jupiter population, in particular type II mi- 1998), or intrusion by the planet into a central cavity or surface 

gration, where the giant planet has grown massive enough to density transition in the disk, decoupling it from the evolution 

open a gap in its protoplanetary disk and migrates inward in of the gas (Lin et al. 1996; Kuchner & Lecar 2002; Masset 

step with the disk's viscous evolution (e.g. Lin & Papaloizou et al. 2006; Papaloizou 2007). Halting migration further out, 

1986; Lin et al. 1996; Ward 1997; Nelson et al. 2000). Giant ex- beyond the < 0.1 AU hot- Jupiter region, may require that gi- 

oplanets at intermediate distances, where eccentricities can be ant planets form late in the lifetime of the gas disk and hence 

high, might be explained by mutual scattering of giant planets only have time for a partial inward migration before stranding 
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at an intermediate distance when the gas is lost (Trilling et al. 
1998). Disks around T Tauri stars are observed to last for ~ 1 
-10 Myr (Haisch et al. 2001) but disperse over a much shorter 
~ 10 5 year timescale (Simon & Prato 1995; Wolk & Walter 
1996): a behavior that may result primarily by accretion of gas 
onto the central star combined with photoevaporative gas loss 
driven by the stellar UV output (Clarke et al. 2001; Alexander 
et al. 2006). Models of this stranding mechanism (Armitage et 
al. 2002; Armitage 2007), which can roughly reproduce the ex- 
oplanet semi-major axis statistics, have raised the possibility 
that fortuitous disk dispersal might also explain the presence of 
the hot-Jupiter population and imply that earlier formed giant 
planets could have been consumed by the cental star. 

If type II migration correctly accounts for the presence 
of the hot-Jupiter population then these giant planets must 
have traversed their inner systems at a time when gas was 
still present and before the completion of terrestrial planet for- 
mation. This prompted initial speculations that such systems 
would be likely to lack any terrestrial planets within their inner 
few AU (Armitage 2003), and since hot-Jupiters are not un- 
common, they have been used to infer significant constraints on 
the abundance of habitable planets (Ward & Brownlee 2000), 
and even their galactic location (Lineweaver 2001; Lineweaver 
et al. 2004). This view is contradicted however by recent mod- 
els that have simulated the process of a giant planet migrating 
through an interior protoplanetary disk (Fogg & Nelson 2005, 
2006; Raymond et al. 2006; Fogg & Nelson 2007a; Mandell et 
al. 2007; Fogg & Nelson 2007b). These find that solid material 
is not predominantly accreted by the giant planet or the cen- 
tral star; instead, solid bodies captured at interior mean motion 
resonances with the giant are shepherded inward an arbitrary 
distance before being randomly scattered into an external or- 
bit. The net result after the migration is a partitioning of most 
of the original disk material into two remnants: a compacted 
remnant interior to the final orbit of the giant, which typically 
accretes in a short timescale to form hot-Earth or hot-Neptune 
planets; and an external disk of scattered bodies. The relative 
predominance of these outcomes has been shown to be sensi- 
tive to the strength of dissipative forces operating at the time of 
migration (Fogg & Nelson 2005, 2007a,b) with scattering be- 
coming increasingly prevalent in late migration scenarios when 
less gas is present. All these studies concur that a scattered disk 
of sufficient mass to support renewed planet formation is likely 
to be generated under a variety of conditions and that terrestrial 
planets should be commonplace in hot-Jupiter systems, rather 
than rare or absent. 

One simplification common to these previous models is that 
the physical mechanism that actually halts giant planet migra- 
tion is not specified or modeled. Type II migration is artificially 
halted when the giant planet has reached a preset final orbit and 
hence is not determined by the structure or evolution of the gas 
disk. Since these models stop migration close to the central star, 
whilst significant gas is still present, they appear most realistic 
in the context of a central gaseous cavity halting mechanism. 
The examples that come closest to implicitly assuming fortu- 
itous gas disk dispersal as the halting mechanism are the late 
scenarios of Fogg & Nelson (2007 a,b) where gas densities have 
fallen to low levels and migration is decelerating (see Fig. 3 in 



Fogg & Nelson (2007a) and Fig. 3 of this paper). However, the 
final hot-Jupiter orbits in these papers are still artificially im- 
posed at 0. 1 AU and are not controlled self-consistently by the 
evolution of the gas. 

Our previous model adopted a 1-D, viscously evolving, gas 
disk algorithm which simulates accretion onto the central star, 
annular gap formation in the vicinity of a giant planet, and self- 
consistent type II migration; for the nebular parameters cho- 
sen, the mass of our gas disk exponentially declined with an 
e-folding time of 582000 years (Fogg & Nelson 2007a). This 
sort of model runs into trouble when simulating the late stages 
of gas disk dispersal as it does not reproduce a final and abrupt 
~ 10 5 year decline that would accord with observations. We 
have corrected this deficiency here by including a photoevapo- 
ration algorithm in our code that gradually erodes and removes 
mass from our gas disk. As shown by Clarke et al. (2001) & 
Alexander et al. (2006), this process has little effect on the evo- 
lution and structure of the gas disk at early times, but comes to 
dominate at later times once the rate of gas loss onto the central 
star due to viscous evolution falls below the photoevaporation 
rate. A rapid dispersal of the remaining gas follows, along with 
the cessation of any ongoing giant planet migration. 

In this paper, we report on the results of a set of scenar- 
ios where giant planet stranding distances are no longer pre- 
scribed but which happen when migration runs out of steam at 
the time of the disappearance of the nebular gas. We therefore 
specifically assume and self-consistently model fortuitous gas 
disk dispersal as the mechanism that finalizes giant planets in 
their post-migration orbits. Terrestrial planetary formation in 
this context is of interest because, for a hot-Jupiter to strand 
at ~ 0.1 AU, it must form and migrate late in the lifetime of 
the gas disk, when gas densities are lower and accretion in the 
inner system is at a more advanced stage than previously con- 
sidered. In addition, a succession of later scenarios than this re- 
sults in a succession of shorter migrations and larger stranding 
distances. This has allowed us to extend the scope of our study 
to model terrestrial planet growth in those 'warm- Jupiter' sys- 
tems that may have originated as the result of a late, partial, 
inward migration. In this paper we define a 'warm-Jupiter' to 
be one orbiting with semi-major axis in the range 0.1 < a < 2.7 
AU, where the outer limit coincides with the snowline. 

The plan of the paper is as follows. In Section 2 we outline 
the additions to our model and the initial conditions of the sim- 
ulations; in Section 3 the results are presented and discussed; in 
Section 4 we consider some caveats, and in Section 5 we offer 
our conclusions. 

2. Description of the model. 

We model our systems using an enhanced version of the 
Mercury 6 hybrid-symplectic integrator (Chambers 1999), run 
as an N + N' body simulation, where there are N protoplanets 
embedded in a swarm of N' "super-planetesimals" - tracer par- 
ticles with masses a tenth of the initial masses of protoplanets 
that act as an idealized ensemble of a much larger number of 
real planetesimals and are capable of exerting dynamical fric- 
tion on larger bodies (e.g. Thommes et al. 2003). The central 
star, giant planet, and protoplanets interact gravitationally and 
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can accrete and merge inelastically with all other bodies. Super- 
planetesimals however are non-self-interacting but subject to a 
drag force from their motion relative to the nebular gas that is 
equivalent to the gas drag that would be experienced by a sin- 
gle 10 km radius planetesimal. Details of these aspects of our 
model are given in Fogg & Nelson (2005). 

We calculate the evolution of the nebular gas using a 1-D 
viscous disk model that solves numerically a modified viscous 
gas disk diffusion equation that includes the tidal torques ex- 
erted by an embedded giant planet (Lin & Papaloizou 1986; 
Takeuchi et al. 1996) and have described its implementation in 
Fogg & Nelson (2007a). The gas responds by depleting over 
time via viscous accretion onto the central star; opening up an 
annular gap centered on the giant planet's orbit; and forming 
a partial inner cavity due to dissipation of propagating spiral 
waves excited by the giant planet. The back reaction of these 
effects on the giant planet is resolved as torques which self- 
consistently drive type II migration. We model the possible ef- 
fects of type I migration (Ward 1997; Papaloizou & Larwood 
2000; Tanaka et al. 2002; Tanaka & Ward 2004; Cresswell & 
Nelson 2006), where a tidal interaction with the gas disk is 
thought to exert an inward radial drift and strong eccentricity 
and inclination damping on protoplanets of ~ 0.1 - 100 M ffl , 
using a simple algorithm described in Fogg & Nelson (2007b). 

As well as modeling the dynamics of gaseous volatiles 
in our model, we also track the movement of presumed solid 
volatiles, such as water ice and hydrated minerals, by labeling 
all particles with a composition based on their original location 
in the disk and summing the composition of protoplanets as 
they grow. We assume a crude three-phase initial radial com- 
position with rocky material originating at < 2 AU, material 
similar to chondritic meteorites between 2 - 2.7 AU, and trans- 
snowline material at > 2.7 AU. We do not assign an actual water 
mass fraction to these phases. 

2.1. Photoevaporation-driven disk dispersal. 

A gas disk that viscously drains onto its central star undergoes 
a power law decline with a long drawn out dispersal in conflict 
with observations that final dispersal occurs over a timescale 
that is short compared with the disk age. Clarke et al. (2001) 
showed that this could be explained by including a model of 
photoevaporation of the disk driven by the diffuse U V flux from 
the central star (Hollenbach et al. 1994). Their results show 
that once the accretion rate onto the star declines to roughly 
equal the outer disk photoevaporation rate, the inner disk ceases 
to be resupplied from larger radii and rapidly drains onto the 
star. The formation of this inner cavity then permits direct UV 
illumination of the outer disk which disperses in turn in ~ 10 5 
years (Alexander et al. 2006). 

For our purposes we need only adopt a simple parameteri- 
zation of this type of photoevaporation model and subtract from 
the right hand side of our disk diffusion equation (Eq. 7 in Fogg 
& Nelson (2007a)) an extra term representing a disk wind: 
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Fig. 1. Upper panel. Mass of the nebular gas vs. time: the blue curve 
represents our former model where mass loss occurs solely via vis- 
cous accretion onto the central star; the red curve represents our new 
model including photoevaporation. Lower panel. Accretion rate onto 
the central star Am, I At. The dotted horizonal line represents the pho- 
toevaporation rate Am w /At = 10" 9 M Q yr~' 



where S w is the rate of change of gas surface density due to 
photoevaporation, r is radial distance, r out is the disk radius, 
and r g is the gravitational radius: the distance beyond which 
ionized gas can become unbound from the star. The constant of 
proportionality K depends directly on the disk's total photoe- 
vaporative mass loss rate m w : 



K = m. 



-2nr 



2.5 



f 

•Jr. 



-1.5 



dr 



(2) 



To fit with our pre-existing viscous gas disk model, we take 
'"out = 33 AU, and adopt r g = 5 AU and m w = 10~ 9 M yr~' 
which gives K = 1.4685 x 10~ 12 g ctrT 2 s _1 . 

In Fogg & Nelson (2007 a,b) we assumed an initial condi- 
tion of a minimum mass solar nebula model (Hayashi 1981), 
scaled up in mass by a factor of three (3xMMSN), extending 
between 0.025 - 33 AU from a solar mass protostar, with an 



initial surface density profile of E g 



and a total mass of 



r g <r < r out , 



0.039 M . Having chosen an alpha viscosity of a — 2 x 10~ 3 , 
we found that, after a short lived ~ 10 5 year period where Z e 
close to the star relaxes to a shallower profile, the mass of the 
gas disk declines predictably with an e-folding time of 582 000 
years. This behavior is illustrated as the upper blue curve in 
the top panel of Fig. 1, which plots the nebular mass vs. time, 
and is compared with the red curve which illustrates the effect 
of including photoevaporation. It is evident that the two mod- 
els only diverge slowly for the first ~ 2 Myr, but thereafter the 
mass of the photoevaporating disk drops steeply and vanishes 
in just a few xlO 5 years. The lower panel of Fig. 1, which plots 
the accretion rate onto the star m„(f), shows that this transition 
in behavior occurs around the time when ra» a; m w . 

The evolution of the gas disk surface density £ e (r, f) is 
shown in Fig. 2, where the uppermost curve represents the ini- 
(1) tial Eg oc r~' 5 profile and the lower curve represent succes- 
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r(AU) 

Fig. 2. Model gas disk surface density evolution. The uppermost curve 
is the initial condition; successive curves are labeled with their age in 
Myr. 

sively more evolved configurations. It can be seen that the evo- 
lution of the nebula speeds up after ~ 2 Myr, with a gap at 
r ~ r g starting to open up at ~ 2.38 Myr, followed by an ac- 
celerated decline of the inner disk thereafter. This behavior is 
qualitatively similar to that described by Clarke et al. (2001), 
with the exception that at late times our outer disk, which is 
truncated to a much smaller radius, is lost even more rapidly. 
This makes no difference to our mechanism for stranding gi- 
ant planets as the divergence occurs when the nebular mass has 
already fallen below the level where it can drive migration. If 
giant planets strand well inward of r g we found that they can 
act against the efficient draining of the last dregs of the inner 
disk onto the central star, slightly altering the picture given in 
Fig. 2. Again however, this effect is minor as it occurs at times 
when the gas is very thin and does not significantly delay the 
date of overall disk dispersal. 

2.2. Initial conditions and running of the simulations. 

We do not consider the earliest stages of planetesimal forma- 
tion and runaway growth in our modeled systems and set the 
t = start date for our simulations to be 0.5 Myr after the start 
of star formation. By this time, we assume that there is no more 
infall of gas onto the protoplanetary disk and the solids compo- 
nent of the inner disk has reached its oligarchic growth stage, 
where a succession of protoplanets, each being a few percent 
of an Earth-mass, have emerged from the planetesimal swarm 
in near-circular orbits and with roughly equidistant spacing in 
units of mutual Hill radii (Kokubo & Ida 2000). We assume a 
central star of 1.0 M with initial surface density profiles for 
gas and solids (Z oc r~ L5 ) taken from a minimum mass solar 
nebular model (Hayashi 1981), which is scaled up in mass by 
a factor of three (3xMMSN) to provide enough mass beyond 
the nebular snowline for a giant planet to form before the loss 
of the nebular gas (Lissauer 1987; Thommes et al. 2003). We 



Table 1. Data describing initial solids disk set-up 





Rocky Zone 


Icy Zone 


Total 




0.4-2.7 AU 


2.7-5.0 AU 


0.4-5.0 AU 


M solld 


9.99 M ffi 


24.65 Me 


34.64 Me 


^proto 


0.025 M e 


0.1 Me 




N 


66 


15 


81 




0.0025 Me 


0.01 Me 




N' 


3336 


2315 


5651 


/proto 


0.17 


0.06 


0.09 



model the gas component of the protoplanetary disk between 
0.025 - 33 AU, giving an initial mass of 0.039 M Q , and com- 
pute its evolution with our photoevaporating viscous gas disk 
algorithm, as described in the previous section and in Fogg & 
Nelson (2007a). We take the alpha viscosity of the gas to be 
a - 2x 10~ 3 , which gives a viscous evolution time at 5 AU ^ 
120000 years. 

We model the solids component of the disk initially be- 
tween 0.4 - 5.0 AU and assume a snowline at 2.7 AU be- 
yond which the mass of solids is boosted by a factor of 4.2 
by the condensation of ices. We generate the initial N-body 
components of the solids disk in an identical manner to that 
detailed in Fogg & Nelson (2005) by starting with initial pro- 
toplanetary masses of 0.025 and 0.1 M ffl interior and exterior 
to the snowline respectively, spaced approximately 8 mutual 
Hill radii apart, with the remainder of the material inventory 
consisting of super-planetesimals with a fixed mass of 10% of 
the initial masses of the local protoplanets. Relevant data for 
the initial solids components are shown in Table 1 which gives, 
for zones interior and exterior to the snowline, values for the 
total mass of solid material M so i;d, the number and mass of 
protoplanets N and m proto , and the number and mass of super- 
planetesimals N' and m s _ p i. The parameter / proto , at the foot 
of Table 1, is the mass fraction of the solids disk contained in 
protoplanets and we use this here as a rough measure of the 
evolution of the disk, taking / pro to = 0.5 to denote the transi- 
tion between oligarchic and chaotic, or 'giant impact', growth 
regimes (Goldreich et al. 2004). 

Our previous approach was to run our combined N-body 
and gas disk model, in the absence of a giant planet, from t = 
to a set of durations distributed between 0.1 - 1.5 Myr, in order 
to mature the disk to different ages and to generate a set of mi- 
gration scenarios. At the end of each of these maturation runs, 
we then introduced a 0.5 Mj giant planet at 5 AU, after remov- 
ing 0.4 Mj of gas from a local disk annulus to provide for the 
giant planet's envelope. This gas is assumed to have accreted 
on top of a 0.1 Mj solid core, composed of material deriving 
from beyond the outer boundary of our modelled solids disk. 
Re-starting the run at this point resulted in an inward type II 
migration of the giant planet which was allowed to continue 
until it reached our prescribed stranding radius of 0.1 AU. This 
approach is not appropriate here as we are specifically assum- 
ing gas disk dispersal as the stranding mechanism and therefore 
need to mature the disk to between the boundaries of a tempo- 
ral window within which our giant planet can both accrete suf- 
ficient gas for its envelope and cease migration at > 0. 1 AU. We 
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Fig. 3. Migration and stranding of 0.5 Mj giant planets in our present 
and previous models. Giant planet semi major axis in AU is plotted 
against t s , time normalized to the start of migration, in units of 10 5 yr. 
Black curves show the behavior of giant planets stranding naturally 
due to gas disk photoevaporation. Numeric labels are launch dates: 
the respective disk ages f/Myr at which the giant planet is introduced. 
Dashed curves show giant planet migrations in Scenarios III, IV, and V 
from Fogg & Nelson (2007a), where migration was artificially halted 
at 0.1 AU. 



located the lower limit of this window, the age where a 0.5 Mj 
giant planet introduced at 5 AU will naturally cease migrating 
and come to rest at ~ 0.1 AU, via experiments with the model 
and found it to be t « 1 .77 Myr. By this time, the mass of the 
nebula has fallen to M gas « 1.25 Mj which is shown as the 
upper dotted line in Fig. 1. The maximum possible upper age 
limit of the stranding window would be when M gas = 0.4 Mj, 
at t as 2.13 Myr (see the lower dotted line in Fig. 1), giving 
a window duration of ~ 17% of the simulated disk lifetime. 
However, this limit would require the unrealistic condition of 
all the remaining nebular gas being accreted by the giant planet. 
Since we do not simulate the process of gas accretion onto the 
giant planet's core, we have arbitrarily restricted the upper age 
limit of the stranding window to t = 1.85 Myr (see the mid- 
dle dotted line in Fig. 1) by which time the mass of the nebula 
has fallen to M gas = 1.01 Mj, reducing the window duration to 
~ 4% of the simulated disk lifetime. We note that the model 
of Armitage et al. (2002) predicts a stranding window of dura- 
tion « 20% of the disk lifetime. This longer duration compared 
to ours is because they adopted a more slowly evolving disk 
model and were simulating the stranding of more massive gi- 
ant planets which migrate more slowly. Clearly, variation of 
the many free parameters in a model such as this can produce a 
variety of stranding behaviors in simulations, but we have not 
considered these here as our main focus is on the effect that the 
giants' stranding distances have on the partitioning of the inner 
system disk and subsequent terrestrial planet formation. 

The behavior of 0.5 Mj giant planets launched into our 
model disks, between the lower and upper age limits discussed 
above, is illustrated by the solid curves in Fig. 3 and shows that 



stranding takes place between 0.13-1.71 AU. Also illustrated 
by the dashed curves in Fig. 3 are the migration trajectories of 
the giant planets in Scenarios III, IV, and V of Fogg & Nelson 
(2007a) where migration takes place in a younger disk and was 
artificially halted at 0.1 AU by a presumed inner disk cavity. 
The fastest of these (launched at 0.5 Myr) takes place at a time 
when the gas disk is still quite massive and completes its mi- 
gration in the viscous evolution time of ~ 120000 years. Later 
scenarios entail longer migration times as the gas mass pro- 
gressively declines and is less effective at driving migration. In 
order for giant planets to strand naturally, our present models 
require still later launch times, in a context where photoevap- 
oration is starting to have a significant influence on the disk. 
Fig. 3 shows that migration speeds are considerably slower and 
decelerate steadily until migration ceases after t s = 600000 
years for the farthest travelling planet and 470 000 years for the 
planet that migrates the least. Migration in all these cases halts 
at t « 2.4 Myr. 1 

Thus, we generated the scenarios for this paper by run- 
ning the model from its initial condition, without a giant planet 
present, to mature the protoplanetary disk to a minimum of 
t = 1.77 Myr and then in successive 10000 year increments 
to t — 1.85 Myr. Given that the reality of strong type I mi- 
gration is controversial, and to bracket the range of possibil- 
ities, two parallel sets of scenarios are generated: one with 
no type I migration forces (Run Set A) and the other with 
type I migration and eccentricity and inclination damping set 
at the maximum rate determined by Eqs. 1 & 2 in Fogg & 
Nelson (2007b) (Run Set B). To accommodate this total of 18 
simulations, a change from our previous scenario ID system 
is also required. Thus, Roman numerals are substituted with 
Arabic numerals and labelling is generated by the following 
formula: Scenario ID = 1 + (f/Myr - 1.77)/0.01. This gives 
Scenarios 1, 2, .., 9 f or t = 1.77, 1.78, .., 1.85 Myr. Run set B 
which includes type I migration is denoted by an I subscript: 
e.g Scenarios li, 2i, .., 9t. During these maturation runs the 
simulation inner edge was set at 0. 1 AU and any material pass- 
ing interior to this boundary was eliminated and assumed to be 
consumed by the central star. The configuration of these ma- 
tured solids disks at t = 1 .77 Myr, the opening of the stranding 
window, are shown in Fig. 4 and Table 2 gives relevant data 
including values for the remaining total mass M so iid, the maxi- 
mum protoplanetary mass m max , the numbers of surviving pro- 
toplanets and super-planetesimals N & N' , and the protoplanet 
to whole disk mass fraction / pro to- Since the scenarios of the 
present work start much closer together in time than those of 
our previous models (0.01 Myr vs. 0. 15 - 0.5 Myr), the state of 
the solids disks generated for Scenarios 2-9 does not change 
greatly from that of Scenario 1 . 

It can be seen, when comparing with the initial condition 
data in Table 1, that planetary growth has been strong, espe- 
cially where no type I migration is operating. In this case (the 
upper panel in Fig. 4), mergers have reduced protoplanets to 



1 Note that the gap between the t = 1.82 and 1.83 Myr curves is 
because of the substantial depletion of disk gas at these late times, 
dictating the removal of the quantity required for the giant planet's 
envelope from a wider annulus of the disk than previously. 
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Fig. 4. Eccentricity vs semi-major axis for 
matured solids disks at 1.77 Myr for no type 
I migration (top panel) and with type I mi- 
gration (bottom panel). Grey dots are super- 
planetesimals; red and blue circles are pro- 
toplanets originating interior and exterior to 
the snowline respectively. Gas densities are 
read on the right hand axes: the upper blue 
lines show gas density at t = and the lower 
cyan lines are the current densities. 



Table 2. Matured solids disk at the start of the stranding window 



fMyr 


1.77 


1.77 


Type I migration ? 


no 


yes 


Msolid 


33.52 M ffi 


30.42 M e 


"*max 


2.67 M e 


1.62 Me 


N 


27 


19 


N' 


1854 


1809 


,/proto 


0.62 


0.52 



a third of their former number, m max is high (there being a 
2.67 M ffl planet present at 2.01 AU), and / plo to indicates that 
the accretion pattern of the disk has progressed way beyond 
oligarchic growth into the chaotic growth regime. Very little 
mass has been lost interior to 0.1 AU (~ 3%) via dynamical 
spreading and gas drag induced orbital decay of planetesimals. 
With type I migration, there is in play an additional preferential 
damping and inward migration of the most massive protoplan- 
ets (clearly visible in the lower panel in Fig. 4) resulting in the 
loss of ~ 12% of the disk mass beyond the simulation inner 
edge. This loss is mostly in the form of large bodies as can be 
inferred from the lower values of m max , N, and f ploto . It might 
be thought that this loss is quite modest considering our inclu- 
sion of type I migration forces. However, in a rapidly dispersing 
gas disk model such as ours, inward type I migration, which is 
proportional to planetary mass, is limited at early times by the 
small size of protoplanets, and at late times by low gas densi- 
ties (see also McNeil et al. 2005; Daisaka et al. 2006); Fig. 4 
shows that at t — 1 .77 Myr, E g interior to 1 AU has fallen by 
two orders of magnitude. The effect of Type I migration on the 
radial distribution of solids disk mass is shown in Fig. 5 where 
the total solids mass for both models at t — 1 .77 Myr is plot- 
ted in 0.5 AU width bins against radial distance. It can be seen 
that beyond ~ 2 AU the radial mass profile of the disks in the 



10-, 




r(AU) 



Fig. 5. Total solids mass in 0.5 AU width bins at / = 1.77 Myr for 
maturation runs with no type I migration (magenta bars) and with type 
I migration (cyan bars). 



two models remains similar, but interior to 2 AU the more rapid 
pace of protoplanetary growth has resulted, in the type I migra- 
tion case, in an inward displacement of mass caused mainly by 
a fractionation of the most massive bodies from the rest of the 
swarm which are now crowding the inner 0.5 AU of the sys- 
tem. It was shown in Fogg & Nelson (2007b) that hot-Earth 
type planets are more likely to accrete and survive when a gi- 
ant planet migrates through such a solids disk, where previous 
type I migration has caused a radial contraction of the inner 
mass distribution. 
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The matured disks summarized above, aged in 10000 year 
stages from t — 1 .77 Myr, are used as the basis for the type II 
giant planet migration scenarios presented here. In each case, a 
giant planet of 0.5 Mj is inserted at 5 AU after removing 0.4 Mj 
of gas from the disk and the inner boundary of the simulation is 
reset to 0.014 AU = 3 R©, which is approximately the radius of 
a solar mass T-Tauri star (Bertout 1989). The giant planet then 
proceeds to clear an annular gap in the gas and undergoes in- 
ward type II migration. The simulations are halted when the 
giant planet strands as a result of the near complete loss of 
the disk gas. In practise, since our viscous disk algorithm re- 
quires a finite amount of gas in each cell to remain stable, we 
assume that the giant halts inward migration when the migra- 
tion rate falls below - 0.2 cm s _1 s - 4.2 x 10~ 7 AU yr _1 . By 
this time the total gas remaining in the entire modeled disk is 
< 10~ 5 M and any error in stranding radius caused by this 
procedure is only on the order of ~ 10~ 3 AU. The symplectic 
time-step for these runs was set to one tenth the orbital period 
of the innermost object which was achieved by dividing each 
simulation into a set of sequential sub-runs with the time-step 
adjusted appropriately at each restart. Since these new simu- 
lations involved a migration of roughly triple the simulated 
duration of our previous models (see Fig. 3), and since small 
time steps were usually needed during the long drawn out 'end 
game' when the giant planet and the shepherded fraction of the 
solids disk are close to their final positions, these runs took a 
particularly long while to complete, requiring 4-6 months of 
2.8 GHz-CPU time each. 

3. Results of the model. 

3.1. System configurations at the stranding point. 

All scenarios, when run to the point at which the giant planet 
ceases migrating, exhibit a varying mix of the same shepherd- 
ing and scattering effects on the solids disk shown in Fogg & 
Nelson (2005, 2007a,b). It is unnecessary therefore to repeat 
the previous procedure of giving a detailed account of the evo- 
lution of one representative case. The big difference is that the 
scenarios presented here result in a more restricted type II mi- 
gration, with the giant planet stranding between semi-major 
axes of fl g a 0. 1 - 1 .7 AU, in a context where all the damp- 
ing forces that are dependent on the disk gas are close to their 
minimum possible values. 

The end points 2 of Scenarios 1-9, those without type I 
migration (Run Set A), are all illustrated in Fig. 6. Their coun- 
terparts, Scenarios It - 9j, with type I migration operating (Run 
Set B), are illustrated in Fig. 7. Comparison of the figures 
shows the familiar partitioning of the solids disk into shep- 
herded interior and scattered exterior fractions, systematically 
truncated by the extent of the traverse of the giant planet. Early 
scenarios (1-3 and h - 3i) are the closest to previous mod- 
els and result in relatively better populated exterior disks and 
sparser interior disks without type I migration in force, and the 
opposite tendency with type I migration. This is in accord with 

2 The phrase 'end point' here refers to the time at which type II 
migration of the giant planet ceases due to gas disk dispersal. 



the previous findings of Fogg & Nelson (2007a) and Fogg & 
Nelson (2007b) respectively. 

No surviving protoplanets are found in or close to the sys- 
tem's maximum greenhouse habitable zone (~ 0.84 - 1.67 AU; 
Kasting et al. 1993) when the giant strands between 0.4 < a g < 
1.2 AU. When stranding occurs at a g < 0.4 AU, late scattered 
protoplanets can find themselves emplaced in the exterior disk 
at distances of < 2 AU and are hence candidates for evolving 
into future habitable planets. However, this eventuality appears 
less likely if type I migration is influential on the dynamics, 
as interior disk fractions evolve closer to the star, are better 
damped, and hence are less likely to lose their contents via 
late scattering. When stranding occurs at a g > 1.2 AU, pro- 
toplanets are found to survive in the inner regions of the HZ 
at ~ 1 AU in both scenario sets. The reason that interior HZ 
planets are found much closer to the giant planet than those in 
external orbits is simply a reflection of the asymmetry between 
shepherding and scattering behaviours. Shepherding of the in- 
terior population occurs between the 2:1 and 4:3 resonances, 
causing material to accumulate via disk compaction between 
a a 0.63-0.83 a g . In contrast, scattering typically results in the 
expulsion of a protoplanet into the exterior disk with an initial 
e > 0.5 and periastron « a g at the time of scattering; hence, ex- 
terior HZ planets are usually found with semi-major axes much 
larger than the final semi-major axis of the giant. Habitable 
planet candidates can therefore be expected if a migrating gi- 
ant planet makes a limited excursion into the HZ. However, if it 
traverses the HZ, such candidates are only expected to be com- 
mon if the giant continues its migration to a radial distance of 
less than half that of the inner edge of the HZ. Whether the po- 
tential habitable planets visible in Figs. 6 & 7 can survive the 
long final phase of accretion that remains to be played out in 
their respective systems is examined in Sect. 3.2. 

The interior disks that result after the giant planet's migra- 
tion stalls show a considerable difference when the two figures 
are compared. When type I migration operates (Fig. 7), inte- 
rior partitions tend to be more massive and have cleared al- 
most all of their planetesimal population. Hence / pr oto ~ 1, but 
large numbers of protoplanets are also found, as type I eccen- 
tricity damping exerted by the residual gas acts to reduce the 
effects of mutual scattering, causing protoplanets to dynami- 
cally settle into stable resonant convoys with closely spaced, 
and near-circular, orbits (originally described by McNeil et al. 
(2005)). For example, the inner eight protoplanets at the end 
points of Scenarios 5i - 7i, from the inside out, are locked 
into a 4:3, 5:4, 4:3, 3:2, 5:4, 7:6, 7:6 configuration of mean 
motion resonances. Oligarchic growth therefore ends in these 
interior systems when the planetesimal field is accreted, but gi- 
ant impact growth is delayed, at least for as long as the disk 
gas persists. When no type I migration operates (Fig. 6), in- 
terior partitions are relatively depleted of mass and contain 
fewer protoplanets, in more excited orbits, alongside a surviv- 
ing population of planetesimals (/ pro to < !)■ In Scenarios 1 - 
2, where the giant planet has experienced the lengthiest mi- 
gration to flg < 0.2 AU, late scattering or accretion has re- 
moved all interior protoplanets (similar to the results of Fogg 
& Nelson (2007a)). In Scenarios 3-5, just one relatively low 
mass hot-Earth remains at the end point. The inner disk fares 
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Stranding point results : with no type I migration 
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Fig. 6. Run Set A. End points of scenarios that exclude type I migration, when the giant planet strands at its final semi-major axis. Eccentricity 
is plotted vs. semi-major axis with symbols colour coded as in Fig 4 and sized according to the mass key. Scenario ID is given at the top left of 
each panel. Protoplanets interior to the giant, or within 1-2 AU, are labelled with their mass in M e . The dotted box shows the habitable zone. 
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Stranding point results : with type I migration 
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Fig. 7. Run Set B. End points of scenarios that include type I migration, when the giant planet strands at its final semi-major axis. Eccentricity 
is plotted vs. semi-major axis with symbols colour coded as in Fig 4 and sized according to the mass key. Scenario ID is given at the top left of 
each panel. Protoplanets interior to the giant, or within 1-2 AU, are labelled with their mass in M e . The dotted box shows the habitable zone. 
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better in late scenarios (7 - 9) when migration is limited to 
a„ > 1 AU. Numerous interior protoplanets remain in these 
cases, and since type I migration forces are absent, accretion 
via giant impacts is not suppressed and protoplanetary growth 
is more advanced. Since both run sets produce such different 
interior partitions, it is of considerable interest to determine if 
this effects their final architectures after the subsequent phase 
of gas-free accretion. This issue is followed up in Sect. 3.2. 

Some features of the interior systems illustrated in 
Figs 6 & 7 are worthy of further comment. 

Scenarios 3-5: Only one interior planet survives in each 
of these scenarios where the giant planet strands between 
fl g = 0.27 - 0.55 AU, and only one of these is found at a 
first order resonance with the giant at the end point. This is 
the 1.90 M ffl planet in Scenario 4, which is captured at the 2:1 
resonance. In the case of Scenario 3, the interior planet is found 
closer to the giant planet, whilst in Scenario 5 the sweeping 2: 1 
resonance has not quite reached the surviving 1 .60 M e planet. 

Scenario 6: The giant planet has stranded at a g = 0.68 AU, 
leaving four surviving protoplanets in the interior partition. The 
outermost of these has an orbit with e a 0.3, close to a 7:4 pe- 
riod ratio with the giant planet, that crosses the orbit of its near- 
est neighbor. To an accuracy of < 1 %, period ratios between 
the protoplanets, from the outside in, are 4:3, 5:4, and 3:2 re- 
spectively, which are a feature reminiscent of the resonant con- 
voys of protoplanets commonly observed in simulations where 
type I eccentricity damping is included. In this case, dynam- 
ical friction from surviving planetesimals exerts the damping, 
but its relative weakness produces an arrangement that is more 
dynamically excited and clearly unstable. 

Scenario 8: This case stands out from its adjacent Scenarios 
7 & 9, where the giant planet also strands at a g > 1 AU, because 
planetary growth in the interior partition appears to be much 
more advanced, resulting in three planets in well-spaced or- 
bits. This is entirely due to chance giant impacts shortly before 
the scenario end point. It is shown in Sect. 3.2 that, if carried 
through into the gas-free phase, accretion within the interior 
partitions of Scenarios 7 & 9 rapidly catches up, with excess 
protoplanets being eliminated by co-accretion or impact onto 
the giant planet. 

Scenario It: The single interior planet resulting from this 
run is the best hot-Neptune/super-Earth analogue generated in 
this paper. Its mass of 7.57 M e puts it within the observed 
range for such objects. However, in contrast with the results 
of Fogg & Nelson (2007b), the planet is not found at the first 
order 3:2 or 2: 1 resonances, but it is located further away from 
the giant planet, at the second order 3:1 resonance (confirmed 
with a plot of librating resonant angles). Study of this system's 
evolution however reveals that this object was originally shep- 
herded inward at the 2:1, but within the last 0.12 Myr of the 
run it scatters with and eventually accretes two other proto- 
planets interior to it, contracting the orbit of the merged body, 
whereupon it is fortuitously captured at the 3:1 resonance. This 
phenomenon of resonance capture through scattering has been 
recently described by Raymond et al. (2008). The arrangement 
of hot-Neptune and hot-Jupiter in Scenario li is most similar 
to that observed between the two inner planets e and b in the 
55 Cancri system (McArthur et al. 2004), although the orbits of 



these two natural objects are separated by a wider ~ 5: 1 period 
ratio. 

Scenario 2\. Two interior planets survive in this case where 
the giant strands at a g = 0.22 AU. The outermost planet of this 
pair is in the 2: 1 resonance with the giant planet. 

Scenario 3t: Two equal mass hot-Earths survive at the end 
point of this scenario (a g = 0.33 AU), but neither of them are 
in resonant orbits. A resonant convoy of five interior planets 
becomes unstable in the last 40 000 yr of the simulation and is 
broken up by giant impacts. 

Scenario 9\. A noteworthy feature of the end point of this 
run is the presence of a 0.27 M e Trojan protoplanet, located 
in a stable co-orbital resonance with the giant planet. This ob- 
ject was captured into this 1 : 1 resonance shortly after the intro- 
duction of the giant and escorts the larger body inward during 
type II migration. Although the capture of the Trojan planet 
in this case may have been influenced by the abrupt introduc- 
tion of the giant planet at the scenario start time, its presence 
at the end point is not necessarily unrealistic or unexpected. 
Dvorak et al. (2004) & Erdi & Sandor (2004) have shown that 
co-orbital motions can be stable for long time scales and 1:1 
capture and co-migration, in a context where type I migration 
forces are active, has previously been observed in simulations 
by Cresswell & Nelson (2006, 2008). 

Alternate visualizations of the data for the end points of 
Scenarios 1-9 and Scenarios 1t-9t are presented in Figs. 8 & 9 
respectively. These data are the equivalent of the 'fate of the 
solids disk mass' results of previous models set out in Tables 
in the previous papers in this series and have to be recast into 
graphical format here since 18 scenarios are under considera- 
tion, rather than 5 or 6. The nine scenarios of each run set are 
indicated on the x-axis with data from Run Sets A & B being 
shown in the left and right hand panels respectively. In Fig. 8, 
the cumulative percentage of the initial solids disk remaining 
is read off the y-axis, with the red bars representing the frac- 
tion remaining in the shepherded remnant, the blue bars that 
in the scattered remnant, and the grey bars that which is lost 

- predominantly via accretion onto the giant planet. In Fig. 9, 
the y-axis gives both / proto for the interior or exterior remnant 
(the symbols being given in the key), and the stranding posi- 
tion of the giant planet a g /AU which is shown by the orange 
line. Horizontal lines f pmt0 = 0.5 and f pmt0 = 1 are drawn for 
reference. 

Fig. 8 shows that the great majority of the solids disk (~ 80 

- 90 %) survives giant planet migration, regardless of whether 
type I migration operates or not, and regardless of the strand- 
ing position of the giant planet. The bias toward scattering be- 
haviour in the absence of type I migration shows well in the 
data for Run Set A as, when a g < 1 AU, less than 10 % of 
the disk mass remains in the shepherded partition. The biassed 
scattering of protoplanets is indicated by the value of exterior 
/proto > 0.5. In contrast, in the presence of type I migration 
(Run Set B), partitioning of the solids disk is largely insensi- 
tive to a g with ~ 30% and ~ 60% of the mass remaining in 
the interior and exterior partitions respectively in all scenarios. 
The bias against the scattering of protoplanets is indicated by 
the value of exterior / proto a; 0.3. 
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Fig. 8. Fate of the disk mass at the end points of Run Set A & B. Data 
include the cumulative % of initial solids disk in the shepherded or 
scattered remnant, or lost from the disk. 
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Scenario 1 



Fig. 9. Protoplanet mass fraction / proto interior to giant (circles) and 
exterior to giant (squares). Also shown is the stranded giant planet 
semi-major axis a e (stepped line) at the end points of both Run Sets. 



One might reasonably speculate, given the details in Figs. 6 
- 9, that the final systems of terrestrial planets that should 
emerge from these scenarios might be either internally or ex- 
ternally weighted (in terms of both planetary numbers and 
masses) depending on the strength of type I migration forces 
operating during gas-phase formation. If the model reflects re- 
ality, its results may help contribute to the debate over the re- 
ality of type I migration once observational techniques have 
advanced sufficiently to enable a more complete inventory of 
terrestrial planets in exoplanetary systems where type II giant 
planet migration is thought to have occurred. 

3.2. Post-migration terrestrial planet formation. 

The results of the model, shown in Figs. 6 & 7, have taken 
the accretion process only to the point at which the giant planet 
strands due to the loss of the nebular gas. Planet formation in all 
exterior and most interior partitions of the original inner system 
disk is clearly incomplete and awaits a much lengthier phase of 
gas-free accumulation. Simulation of at least part of this gas- 
free accretion phase is of clear interest to see if the conclusions 
drawn from running the model to the end point of migration 
might stand up over the long-term. Would a more complete ac- 



cretion reduce the observed differences between the two run 
sets? What will be the fate of the closely packed resonant con- 
voys of protoplanets created by type I migration (Fig. 7)? Are 
habitable zones containing protoplanets at the end of the mi- 
gration epoch still likely to contain planets over the long-term? 

Thus, after switching off all residual forces exerted by the 
last traces of gas, simulation of all the scenarios of both Run 
Sets has been extended with the aim of running them to a sys- 
tem age of 30 Myr (f = 29.5 Myr). However, this is not prac- 
tical in every case as the time-step is constrained by the orbits 
of the innermost objects. The method adopted here is to choose 
a time step equivalent to one tenth of a circular orbit with a 
semi-major axis equal to the periastron of the innermost large 
object. This strategy typically dictates a time-step of 2 - 4 days 
and serves well in all but two scenarios. These are Scenario 3 
in which a 2 day time-step inadequately resolves the orbits of 
a planetesimal swarm interior to the hot-Earth at 0.18 AU, re- 
sulting in an unrealistic loss of these objects via accretion onto 
both the star and the hot-Earth; and Scenario li where a hot- 
Neptune at 0.065 AU dictates too small a time-step - a prob- 
lem that is tackled by splitting the problem into sub runs - 1) 
the hot-Neptune and hot- Jupiter with a 0.3 day time-step, 2) the 
hot- Jupiter and external disk with a 2 day time-step, and then 
assembling both as a composite. Even so, eight months after 
beginning most of these extended runs it has not been possible 
before submission of this paper to complete every simulation. 
All scenarios in Run Set A have attained the goal of 30 Myr, but 
in Run Set B, where previous type I migration has typically re- 
sulted in protoplanets closer to the star at the point of gas loss, 
the results presented here vary in their ages from 20 - 30 Myr. 

The results of these extended simulations are illustrated 
in Figs. 10 (Run Set A) & 11 (Run Set B) and should be 
compared with their migration end-point equivalents in Figs. 6 
& 7 respectively. A number of general observations can be 
made. 

1. Accretion interior to the giant planet. Solitary, or paired, 
hot-Neptunes or hot-Earths that are present at the end of the 
migration epoch (e.g. Scenarios 3-5 and Scenarios li - 3i) are 
found to survive the extended gas-free accretion phase with lit- 
tle change in their orbital parameters. The crowded interior disk 
partitions that result in later scenarios where the giant comes 
to rest at greater distances from the central star undergo rapid 
evolution to a state of near completion, with giant impacts thin- 
ning down the number of protoplanets which mop up almost all 
residual planetesimal debris. 

As discussed in Sect. 3.1, giant impact growth in Run Set 
A is already underway before the end of the migration phase 
as the protoplanets within interior partitions are sufficiently dy- 
namically excited to exhibit crossing orbits. This process com- 
pletes rapidly after gas loss, with protoplanetary mergers re- 
ducing their numbers by ~ 50%, resulting in apparently stable 
multiple interior systems of two or three hot-Earths separated 
by > 20 mutual Hill radii (note Scenarios 6-9, Fig. 10). 

In Run Set B, type I migration forces have previously 
suppressed giant impact growth, causing protoplanets to 
stack into crowded resonant convoys (e.g. Scenarios 5i - 9i, 
Fig. 7). After the gas is gone in these cases, all eccentricity 
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Gas-free accretion results 
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Fig. 10. Run Set A. Gas-free accretion results of scenarios that exclude previous type I migration. Eccentricity is plotted vs. semi-major axis 
with symbols colour coded as in previous examples and sized according to the mass key. Scenario ID is given at the top left of each panel. 
System age is given at the top right of each panel. Protoplanets interior to the giant, or within 1-2 AU, are labelled with their mass in M ffi . The 
dotted box shows the habitable zone. 
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Gas-free accretion results. 
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Fig. 11. Run Set B. Gas-free accretion results of scenarios that include previous type I migration. Scenario ID is given at the top left of each 
panel. System age is given at the top right of each panel. Protoplanets interior to the giant, or within 1-2 AU, are labelled with their mass in 
Mm. The dotted box shows the habitable zone. 
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Fig. 12. Destabilization of the resonant convoy in Scenario 5i. 
Temporal evolution of the periastron, semi-major axis, and apastron 
for each object is shown, with the giant planet drawn in orange at the 
top of the graph. Protoplanets undergo a sequence of closest neighbour 
mergers, reducing their number from eight down to a stable quartet of 
hot-Earths. 
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Fig. 13. Inner systems of Scenarios 7j & 9i run to t = 100 Myr. 

scattered planetesimal population has resulted in a degree 
of circularization of the orbits of scattered protoplanets, 
especially in the case of Run Set B, where dynamical friction 
is particularly strong due to the low exterior value of / pro to (see 
Fig. 9 and associated discussion). 



damping ceases and the convoys start to destabilize within a 
few Myr, typically by protoplanets merging with their nearest 
neighbours, ultimately reducing their numbers by ~ 50%. This 
occurs fastest in convoys that are the most compressed by the 
position of the giant planet (i.e. restricted within the narrowest 
annulus) and the process is clearly illustrated in the case of 
Scenario 5i in Fig. 12. Generally however, systems where 
strong type I migration has operated in the gas phase still retain 
more interior planets, although is it not certain that all giant 
impact growth is completed within the interior partitions of the 
later scenarios illustrated in Fig. 11. Scenarios 7i & 9i have 
interior protoplanets still separated by as little as ~ 10 mutual 
Hill radii and some traces of the original resonant convoy 
structure remains, such as the 4:3, 5:4, 4:3 commensurability 
between the inner four planets of Scenario 7t. To test the 
stability of the still crowded inner systems shown in Fig. 11, 
the orbits of the giant and the remaining interior planets in 
Scenarios 7t & 9i were integrated to t — 100 Myr and the 
results are displayed in Fig 13. It can be seen that whilst the 
inner system of seven planets in Scenario 9i is stable over this 
extended time frame (including the orbit of the Trojan planet), 
that of Scenario 7i undergoes further orbital evolution and 
accretion resulting in a system of four survivors. Whatever the 
final outcome in these cases however, hot and warm-Earths, 
including multiple systems, are predicted as an outcome of 
type II migration induced disk compaction in all but the 
earliest scenarios of Run Set A. 

2. Accretion exterior to the giant planet. Planetary growth in 
exterior partitions at > 2 AU is ongoing and not yet complete. 
This is to be expected as dynamical times are longer, scattered 
material is spread over a larger volume, and does not have 
an opportunity to interact with disk solids originating beyond 
5 AU which have not been modelled. However, it can be 
seen in Figs. 10 & 11 that dynamical friction exerted by the 



3. Planetary occupants of the habitable zone. Three planets 
in each run set are found in their system's habitable zone 
(~ 0.84 - 1.67 AU) after the extended gas-free runs. In Run 
Set A, these are 3.25, 5.04 & 2.00 M e planets in Scenarios 1, 3 
& 9 respectively (Fig. 10); and in Run Set B a 1.66 M® planet 
in Scenario 8i and 1.26 & 0.27 M e planets in Scenario 9i 
(Fig. 11). This latter object is the Trojan planet captured into a 
1 : 1 resonance with the giant early on in the migration phase 
(discussed in Sect. 3.1) and dragged inward to a — 1.66 AU, 
just inside the outer edge of the habitable zone. This exotic 
world remains in a continuously stable orbit up to a system 
age of 30 Myr, and has the potential for a prolonged existence 
(Dvorak et al. 2004; Erdi & Sandor 2004 Fig. 13). The finding 
articulated in the previous Section that the giant planet must 
make either a limited excursion into the habitable zone, or a 
complete traversal down to < 0.5 x the radial distance of the 
inner edge of the HZ, for a habitable planet to be possible 
is reinforced by the results of the extended simulations. 
Potentially habitable planets are found in systems where the 
giant planet lies outside the region of ~ 0.3 -1.5 AU, except 
in the case where a Trojan planet accompanies a giant stranded 
within the HZ. 

4. Volatile endowments. In agreement with previous models 
(Fogg & Nelson 2007 a,b) and the work of Raymond et al. 
(2006) and Mandell et al. (2007), the habitable planets gen- 
erated in the present model are predicted to be richly endowed 
with volatiles driven inward from beyond the nebular snow- 
line by the migrating giant planet. The composition of the six 
habitable planet candidates, using the same crude three phases 
described in Fogg & Nelson (2007a), are shown as pie charts in 
Fig. 14 where the red sectors refer to dry material originating 
interior to 2 AU, the green sectors to chondritic material origi- 
nating between 2 and 2.7 AU, and the blue sectors indicate the 
fractions originating from beyond the snowline at > 2.7 AU. 
Even with no previous type I migration operating (Run Set A), 
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Fig. 14. Material composition of habitable planet candidates from our 
simulations. Red indicates material of provenance a < 2 AU; green 
indicates material of provenance 2 AU < a < 2.7 AU; and blue indi- 
cates trans-snowline material originating from a > 2.7 AU. The top 
row of charts are the results of Run Set A and the bottom charts the 
results of Run Set B. 

the habitable planet candidates mentioned above have incorpo- 
rated roughly a third to a half of their material from beyond the 
snowline. Where there has been pre-existing type I migration 
(Run Set B), more like three quarters of the planets' mass orig- 
inates from beyond the snowline. Whether the habitable planets 
in question lie interior or exterior to the giant seems to make lit- 
tle difference to the composition. All such worlds are predicted 
to be richly endowed with water and gases. 

4. Discussion. 

The key finding of the model presented in this paper, and 
its predecessors (Fogg & Nelson 2005, 2006, 2007a,b), that 
giant planets do not remove the portion of a protoplanetary 
disk they migrate through, or prevent subsequent terrestrial 
planet formation, is robust in the sense that every model 
variant, and independent work (Raymond et al. 2006; Mandell 
et al. 2007), has reproduced it. However, these models have 
only explored a modest region of parameter space relevant to 
the problem, and have inevitably adopted assumptions that 
simplify or omit potentially important physical processes. 
We have discussed a number of issues concerning parameter 
variations and neglected physical processes in our previous 
papers and restrict the discussion here to three caveats not 
previously considered. 

ij Solids disk inner boundary. The inner system solids disk 
adopted here begins with its inner boundary set at 0.4 AU. Its 
contents are then free to evolve inward during the maturation 
runs preceding the introduction of the giant at the start of a 
migration scenario. By this time, in low dissipation models 
(Run Set A), dynamical spreading typically results in a modest 
inward movement of the disk edge to ~ 0.3 AU (see the upper 
panel of Fig. 4) so that when the giant planet subsequently 
migrates down to 0. 1 AU it passes through the entire inner disk 
annulus. This is not necessarily the case in high dissipation 
models (Run Set B) where greater inward movement of the 



solids disk contents occurs during the maturation phase (see 
the lower panel of Fig. 4). Thus, the question is raised as to 
whether the lack of interior hot-Earths resulting in Scenarios 1 
& 2 is a realistic finding or an artifact of the initial position of 
the solids disk inner boundary. What if this boundary had been 
set at 0. 1 AU instead? This would have provided an additional 
~ 3 M ffl of material from which potential interior planets 
could have been assembled. It cannot be ruled out therefore 
that some hot-Earths could grow and survive even in early 
scenarios of Run Set A, especially if the giant planet stops 
short of their position. However, one would still expect them 
to be of lower mass and more prone to orbital destabilization 
than in high dissipation models where eccentricity damping is 
in play and where there has been much more extensive prior 
inward migration of solids. 

ii) Random component of initial conditions. The solids disk 
with which we begin our simulations at t — conforms as a 
whole to a 3xMMSN model, but is unique in the sense that 
its individual particles start with randomized orbital elements. 
This raises the question of how the conclusions of this paper 
might depend on the seed number used to randomize these 
initial conditions. We have not tested the influence of alternate 
seed numbers on the particular set of simulations presented 
here. However, comparison with our previous work, where 
different integers have been used to seed the initial setup, 
shows that in scenarios where disk maturity and dissipative 
forces are roughly equivalent, the overall architectures of the 
post-migration systems are similar even though details of 
individual planets are different. This is because the general 
physical processes shown to characterize the migration of a 
giant planet through an inner system disk, which lead to the 
division of the original disk contents into internal and external 
remnants, are largely independent of the random component 
in the initial conditions and are influenced primarily by the 
level and variety of dissipative forces in play - i.e. by the 
maturity of the solid body swarm, the density of the gas, and 
whether type I migration is active. An example of such a 
comparison can be seen by referring to Fig. 12 from Fogg & 
Nelson (2007a), where the results of a late scenario, matured 
for 1 .5 Myr before the onset of giant planet migration, in the 
absence of type I migration, are plotted. The lower panel of 
this figure shows the evolution of this model 2 Myr after the 
giant planet has stranded at 0.1 AU and its overall system 
architecture is similar to its closest equivalent in the present 
work: Scenario 1, shown in the top panel of Fig. 6. No interior 
hot-Earth is present and protoplanets amounting to several M ffl 
have been externally scattered to between 1-3 AU in both 
cases. The post migration evolution of this hot- Jupiter system 
from Fogg & Nelson (2007a) has been run further to 30 Myr 
and the results illustrated in Fogg (2008). The resemblance 
between that run and the matured Scenario 1 in this paper 
(Fig. 10, top panel) is even closer with a ~ 3.3 M ffi planet being 
present in the habitable zone of both systems. We conclude 
that although the specific properties of the planetary systems 
which arise from our simulations depend on the exact details 
of the initial setup, the general trends that we observe in our 
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outcomes are robust. 

Hi) Ocean planets. One effect of giant planet migration is to 
drive large quantities of icy material into the inner system, with 
the outcome that terrestrial planets that form in the aftermath 
are predicted to be richly endowed with volatiles and enveloped 
by global oceans. There is nothing about 100% ocean cover that 
necessarily rules out the presence of life. Indeed, life may have 
started in the Earth's oceans, and exploration beneath ocean- 
covered Europa's ice shell has long been regarded as a top pri- 
ority by astrobiologists (Reynolds et al. 1983, 1987). Ocean 
planets in the ~ M e range however have a problem not shared 
with Europa: if the depth of their oceans is > 100 km, their 
floors are composed of high pressure phases of ice, rather than 
rock, with the ice-silicate interface (and much potentially cru- 
cial chemistry) sealed off beneath an icy mantle 100s or 1000s 
of km thick, depending on the planet's total mass and bulk wa- 
ter content (Leger et al. 2004). How this would effect the so- 
lute content of the oceans is unknown but one might specu- 
late that some trace elements essential for terrestrial life might 
be lacking 3 . Meteoritic infall after planetary differentiation is 
complete would re-supply some rock-forming elements for dis- 
solution into the ocean, but solid material that deposits on the 
ocean floor would be denser than the ice beneath and prone to 
removal by sinking into the mantle. A related problem deepens 
this uncertainty. Planets over-endowed with water might sim- 
ilarly possess massive atmospheres with unexpected composi- 
tions and properties. If ocean planets are more like thawed out 
versions of Titan, rather than water-rich Earths, then habitable 
zone calculations that rely on models of modest CO2/H2O at- 
mospheres, coupled to a carbonate-silicate cycle and continen- 
tal weathering (e.g. Kasting et al. 1993; Selsis et al. 2007; von 
Bloh et al. 2007), may not be relevant. 

If the water content of trans-snowline planetary building 
blocks is assumed to be ~ 75% (~ 1 — Z" 1 where /; ce = 4.2 
is the factor in the MMSN model by which nebular condensate 
masses are enhanced due to ice condensation), then the mod- 
els presented here predict HZ planets composed of ~ 20% - 
60% water. This is a vast inventory compared with the Earth's 
~ 0.1% water. It is likely though that these figures are over- 
estimated. If the snowline is further from the star than the 
2.7 AU assumed in the MMSN Model, or if the giant planet 
forms closer to it, the material shepherded inward by type II 
migration will have a lower volatile content. In the models of 
Raymond et al. (2006) and Mandell et al. (2007), where the 
snowline is placed at 5 AU, HZ planets are produced with a wa- 
ter content of ~ 10% - a reduction from the above estimate, but 
still a factor of 20 greater than the typical outcome of ~ 0.5% 
water content when models are run without a migrating giant 
planet (Raymond et al. 2004, 2007). However, all these mod- 
els almost certainly overestimate the quantity of water retained 
by growing planets since loss of volatiles during accretion is 
not accounted for. Extensive depletion of both atmosphere and 
ocean could result from giant impacts late in formation (Genda 
& Abe 2005; Asphaug et al. 2006; Canup & Pierazzo 2006) 



3 The elemental requirements of non-terrestrial life are another 
matter and unknown. 



which could strip nascent ocean planets down to a more Earth- 
like remnant. All of the HZ planets generated by the models 
presented here undergo at least one giant impact at some point 
in their evolution, with a minority suffering a high velocity col- 
lision with a comparable sized object during the later stage of 
accretion within the scattered disk. Thus, whilst ocean planets 
in habitable zones are a robust prediction of the model, it is 
premature to say that they would invariably occur in nature and 
would inevitably be uninhabitable. 

5. Conclusions. 

All previous published models of terrestrial planet formation in 
the presence of type II giant planet migration have neglected 
the issue of what causes the giant planet to come to rest in 
its final orbit and have simply switched off migration forces 
at a pre-defined distance from the central star (Fogg & Nelson 
2005, 2007a,b; Raymond et al. 2006; Mandell et al. 2007). This 
latest model is the first simulation of this type to include a self- 
consistent scenario for the stranding of the giant planet, which 
comes to rest naturally at the point when the gaseous fraction 
of the protoplanetary disk is lost via a combination of accretion 
onto the central star and photoevaporation. 

The results of the earliest scenarios - those which simu- 
late hot-Jupiter emplacement at a g « 0.1 AU - are very similar 
to the late scenarios of Fogg & Nelson (2007a,b), where the 
evolution of the nebula is advanced and dissipative forces that 
stem from the presence of gas are relatively weak. Most of the 
inner system disk solids survive the traverse of the giant planet 
and most of this surviving matter is scattered into higher or- 
bits where planet formation can resume. Water-rich habitable 
planets are possible within a habitable zone that is far removed 
from the final position of the giant planet. Where dissipative 
processes are enhanced by type I migration, one or two hot- 
Earths or hot-Neptunes are found to persist in orbits interior to 
the giant planet. 

In 'warm-Jupiter' type scenarios, where a g » 0. 1 AU, the 
bulk of the disk solids survive as before, but more complex sys- 
tems of hot terrestrial planets interior to the giant are predicted, 
in larger numbers with increasing a g . If strong type I migration 
forces are a genuine influence on planet formation, relatively 
crowded systems of interior planets are predicted, although it 
has not been practical to check this forecast with Gyr-long in- 
tegrations. Habitable planets in low eccentricity warm- Jupiter 
systems appear possible if the giant planet makes a limited 
incursion into the outer regions of the habitable zone, or tra- 
verses its entire width and keeps going until it ceases migrat- 
ing at a radial distance of less than half that of the HZ's inner 
edge: i.e. the giant planet's final orbit lies outside the region of 
~ 0.3 - 1.5 AU for a 1 M star. If giant planets strand within 
this region, the presence of habitable planets remains a possi- 
bility if they are located in a stable 1 : 1 resonance with the giant. 
These findings should hold for stars with different luminosities 
with an appropriate scaling of radial distance. 

However, warm-Jupiter type exoplanets with near-circular 
orbits, i.e. e < 0.05, are relatively uncommon and e is typically 
observed to be considerably higher. Many of these exoplanets 
may have had an origin involving mutual giant planet scatter- 
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Fig. 15. Data of low eccentricity exoplanets using data from the 
Extrasolar Planets Encyclopedia: http://exoplanet.eu. Plotted objects 
correspond to the type of warm-Jupiters within the remit of the 
model presented here: solitary giant planets orbiting single stars with 
0.1 AU < a < 1.8 AU, e < 0.05. Orange symbols represent exoplanets 
within main sequence star systems; red symbols represent exoplanets 
orbiting post-main sequence stars. 

ing, perhaps combined with migration (e.g. Ford et al. 2001; 
Papaloizou & Terquem 2001; Moorhead & Adams 2005), as 
opposed to the strongly damped pure type II migration mech- 
anism adopted here. Nevertheless, low-e and solitary warm- 
Jupiters are known, which presumably could have originated 
in a similar manner to that simulated here. Fig. 15 uses data 
from the Extrasolar Planets Encyclopedia to plot exoplanetary 
data between a = 0.1 — 1.8 AU vs. e — - 0.05 which might 
be regarded as lying within the remit of the model. Solitary 
giant planets around single stars are plotted and labelled with 
their names. Four of these exoplanets, p CrB b, HD 27894 b, 
HD 192263 b and HD 4208 b occur in systems where there is 
dynamical room for terrestrial planets in the habitable zone ac- 
cording to the stability calculations of Jones et al. (2005, 2006). 
Systems containing the exoplanets HD 224693 b (Jones et al. 
2006), HD 104985 b and HD 59686 b (Jones et al. 2005) could 
also have hosted planets in their habitable zones in the past be- 
fore their primaries left the main sequence. According to the 
results of the model presented in this paper, it appears feasible 
that habitable planets could have originated in these systems as 
well. 
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